Convective dust clouds in a complex plasma 
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The plasma is generated in a low frequency glow discharge within an elongated glass tube oriented 
vertically. The dust particles added to the plasma are confined above the heater and form counter- 
rotating clouds close to the tube centre. The shape of the clouds and the velocity field of the 
conveying dust particles are determined. The forces acting on the particles are calculated. It is 
shown that convection of the dust is affected by the convective gas motion which is triggered, in 
turn, by thermal creep of the gas along the inhomogeneously heated walls of the tube. 



Convective motion of micro particles in complex 
(dusty) plasmas is a phenomenon that is often observed 
in very different experimental conditions. In particular, 
vortices in complex plasmas can be produced both in 
ground-based laboratories and under microgravity con- 
ditions, in dc and rf discharges of fairly different con- 
figuration, upon inhomoffeneous heating and in rather 
isothermal environment 

There are many publications, both theoretical and ex- 
perimental, in which the origin of the vortex motion in 
complex plasmas has been investigated [2I, [sl, 0, [H, 0, 0, 
[s, 9]. Basically, there are two mechanisms that can pro- 
duce vortices: This is either the presence of non-potential 
force(s) exerted on charged micro particles in the dis- 
charge (due to inhomogeneous charges [7] or because of 
ion drag [9]) or the convective motion of the background 
neutral gas or dust particles themselves [10]. This clearly 
indicates that the nature of such vortices - despite of 
quite similar appearance - might be very different as well. 

In this paper we report on a recent experiment per- 
formed in a low frequency glow discharge under gravity. 
The PK-4 setup ll|,ll2| was used, which represents a typ- 
ical configuration of a dc discharge employed to study 
complex plasmas. We investigate convective motion of 
micro particles in the presence of an external controllable 
heating and unambiguously show that under such condi- 
tions vortices in complex plasmas occur due to the neu- 
tral gas convection. Clouds of micro particles in this case 
resemble convective clouds in the atmosphere produced 
from warm air pockets rising upwards and composed of 
water droplets, ice crystals, ice pellets, etc. 

The experiments were conducted in a complex plasma 
produced by a low frequency discharge (LFD) in Neon 
gas of the PK-4 facility |l2[. The plasma chamber 
consists of an elongated glass tube as shown in figure [H 
To remove all rest gases which might disturb the exper- 
iments the chamber was evacuated for several hours to 
reach a base pressure below 10~^Pa. Neon gas with con- 
stant pressure between 30Pa and lOOPa were used during 
the experiments. 

A gas discharge was maintained between the electrodes 
by a regulated discharge current of 1mA. The discharge 
voltage polarity was changed with a frequency of Ik Hz 
(50% duty cycle). This frequency is more than an or- 
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FIG. 1: Experimental setup. 



der of magnitude higher than the dust-plasma (response) 
frequency [1], so that the effective longitudinal force on 
micro particles (electric plus ion drag) vanishes. 

In course of the experiments micro particles (spherical 
Melamine- Formaldehyde particles with a mass density of 
l.blg/cm^) of two different radii = 1.64/im or 3.05/im 
were injected by a dispenser from above into the plasma 
and confined inside of the vertical glass tube. The mi- 
cro particles were illuminated along the tube axis by a 
laser sheet (wave length 686nm, power 20mVK, thick- 
ness 100/im ± 30/im, width 15mm) and the images were 
recorded by a PCO 1600 camera. 

A metal wire (alloy of 70% Ni, 11% Fe, and 14% Cu 
and 0.5mm in diameter) was put around the lower part 
of the tube to produce a temperature gradient. A dc 
current between 0.5 A and 1.8 A was applied to this wire 
using a lOF power supply. For avoiding unwanted elec- 
tromagnetic effects on the plasma only one loop of the 
wire was attached. The wire was insulated against direct 
contact to the tube walls. The resulting temperature dis- 
tribution was simultaneously measured at eight points on 
the glass tube at and above the heating wire to extract 
the temperature gradient over a range of 8cm. 

Due to the gas temperature gradient a thermophoretic 
force acts on the micro particles \1^. This was used to 
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FIG. 2: Measured temperatures (solid lines) and calculated 
temperature gradients (dashed lines). 1: temperature and 
temperature gradient for Vp = 1.64/im particles, 2: same for 
rp = 3.05/im particles. 

levitate micro particles against gravity in complex plasma 
experiments [1^. The thermophoretic force acting on 
the micro particles is given by 

Fth = -IthVTn, (1) 

where 7t^[N • cm/K] ^ 2.17 x 10"^ (rp[/im])2 [14]. The 
thermophoretic force required to levitate particles is 
Fii^ = rripg 2.74 x 10~^^A^ for the smaller ones and 
0:^ 1.76 X 10~^^A^ for the larger ones {rup is the particle 
mass). Therefore according to ^ temperature gradi- 
ents of \dT/dz\ :^ A.dbK/cm for the smaller particles and 
S.64:K/cm for the larger ones are required to compen- 
sate gravity. 

The measured temperatures along the tube and the 
temperature gradients are shown in figure O where z = 
is the position of the heating wire. It turns out that the 
measured maximum values of the temperature gradients 
are smaller than those necessary for thermophoretic lev- 
itation. Nevertheless we observe that the particle cloud 
is confined in the area above the heating wire. Figure [3] 
shows vertical cross sections through the center of the 
clouds obtained for different pressures. As we pointed 
out above, plasma forces cannot be responsible for this 
effect: because of the fast polarity switching, neither a 
longitudinal electric force nor an ion drag force [15] could 
support the levitation. In addition, an intense convective 
motion of the particles is observed within the levitated 
clouds. 

The micro particles are frictionally coupled to the neu- 
tral gas. Therefore, it is natural to assume that gas con- 
vection is induced by the local tube heating and that 
the micro particles are dragged by the gas motion. Fa- 
vorably directed, the gas flow could help to levitate the 
particles in addition to thermophoresis and also could 
result in global dust rotation. Note that F^^ alone can- 
not cause the convection, because it is a potential force. 



VxFth = 0. 

The neutral drag force acting on a spherical particle 
within a gas flow is given by [l6[ 

Fn =7n(v/ - Vp), (2) 

where the Epstein friction coefficient reads 
7n[N-s/cm] 0^ 2.7 x 10"^^ (rp[/im])2p[Pa] for Neon 
at room temperature, = 300 K. 

The radial electric field of discharge exerts an addi- 
tional force Fr [17] which confines the particle cloud in 
radial direction. We observed that the micro particles 
after switching off the discharge expand across the tube, 
exploiting practically the entire convective field. Par- 
ticles keep rotating one or two cycles before eventually 
falling down. 

As an example, let us extract the velocity field vj of 
the convective gas flow and the radial electric force F^ 
for the experiment with 3.05 jam particles at a pressure 
of 50Pa. 

Particle motion is recorded at a frame rate of 500 fps in 
course of the experiment comprising two separate cases: 
with and without plasma. The particles are detected in 
each frame and then based on the position of each particle 
in a few consecutive frames their velocities and acceler- 
ations are extracted. Based on this data, the particle 
velocity profiles are reconstructed in the entire convec- 
tive dust cloud. Particle velocities and accelerations are 
derived by fitting cubic splines to the complete particle 
trajectories. 

The trajectories of individual particles are determined 
by the balance of forces: 

mpT = Fth + Fn + F^ + rripg. (3) 

To separate the remaining unknown forces - the radial 
electric force and the neutral drag force - appearing in 
(|3]) , we analyze first the data when no plasma was in the 
tube, and, hence, no radial electric force {Fr = 0) acts on 
the particles. This allows us to get the velocity field of the 
conveying gas molecules. Next with this information and 
assuming that switching on and off the plasma does not 
alter the gas convection, we calculate the radial electric 
force of the plasma in the region where particle trajecto- 
ries recorded with and without plasma overlap (cf. figure 
Hj). Switching off the plasma is a very fast process com- 
pared to the time particles spend conveying in the tube. 
Particles trajectories for both cases are shown in figure 
m In figure [Ha) the particle cloud in the presence of the 
plasma is a bit higher and further away from the tube 
wall. After switching off the plasma [Hb) particles come 
closer to the wall and the cloud expands further down. 
The area of particle clouds and the direction of rotation 
for these two cases are represented in (c) with the over- 
lapping region marked. For the overlapping region any 
differences in the particle motion between the cases with 
and without plasma must be due to electrical forces as 



FIG. 3: Dust clouds of rp = 1.64/im particles at 30, 50, and lOOPa (from left to right). The dash-dotted Unes indicate the 
center of the tube. The field of view is 21 x 26 mm. 




FIG. 4: Convective dust clouds of rp = 3.05/im particles at a pressure 50Pa for two different cases, (a) with switching on of the 
plasma and (b) with switching off the plasma, (c) The superimposed cloud positions for these two cases. The arrows indicate 
the directions of particle rotations. The region of overlapping trajectories is shaded. The dashed lines indicate the center of 
the tube. 



the distribution of temperature and free gas convection 
are not influenced by the discharge. 

The neutral drag force on the particles depends on the 
relative velocities of the particles and gas. This is the 
reason why the particles do not follow the gas flow tra- 
jectories but rotate eccentrically from the center of gas 
convection in the area of upward gas draft. Figure [5] 
presents the velocity field of the free gas convection as 
calculated for the case without plasma. The particle tra- 
jectories are overplotted. The distribution of flow veloc- 
ities has a clear rotational tendency. 

The radial force on the particles when plasma is on (cf. 
figure [6]) can now be fitted using the gas flow velocities 
in the overlapping region. The radial force is zero at the 
tube axis and rapidly increases towards the tube walls. 

The results obtained above clearly demonstrate the 
presence of gas convection. This is not free convection, 



though, merely because the onset of free convection is 
too high in terms of the critical Reyleigh number (see, 
for instance [l^) to cause gas flow under the conditions 
of our experiments. It is well known, however, that if one 
puts a non-uniformly heated body in a rarefied gas, the 
gas starts moving along_the body in the direction of the 



temperature gradient p^, [191, |20[ • This phenomenon is 



referred to as thermal gas creep (or thermal gas slip). It 
was predicted theoretically by Maxwell [2l[ and verified 
experimentally by Reynolds [22[ and is governed by the 
relation 



Vtc = Ktc^Vw \nTw, 



(4) 



where Vtc is the velocity of creep at the body sur- 
face, V is the kinematic viscosity of the gas, Tw is the 
temperature of the heated body, i.e. the glass walls in 
our case, and || indicates the component of the gradi- 
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FIG. 5: Averaged gas flow velocity field (vectors) superim- 
posed with particles trajectories. The dashed line correspond 
to the center of the tube. 
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FIG. 6: Reconstructed radial force for the overlapping region. 
The dashed line indicate the center of the tube. 



ent along the surface. For a long tube, the radial dis- 
tribution of the (longitudinal) velocities is well known, 
"^z = VTc(2r^/i^^ — 1) [20]. In our geometry the ver- 
tical temperature gradient is negative dTw/dz < (see 
figure [2]). Hence, the gas should flow downwards along 
the tube walls with Vz{R) = Vrc and upwards near the 
tube axis with ^^(0) = —Vrc, which is in agreement 
with our observations (see figure [5]). For the quantita- 
tive comparison with the experiment, we rewrite Eq. (|4]) 



in the following form: Ktc = I ^tc /^V;^ In Tv^|. Based 
on the experimental data shown in figures [2] and [5] we 
get Ktc — 1, which is in agreement with theoretical ex- 
pectations. (The theoretically possible range of values of 
Ktc amounts to 0.7-1.2 [l9^). Note also that the magni- 
tude of the convection velocity decreased monotonically 
with pressure, which is also in line with the theoretical 
estimates: combining Eqs. (|3]) and (|4]), we obtain that 
the velocity should scale as (xp~^. 

To conclude, we observed levitation of a particle cloud 
in a vertical glass tube, above a heated wire. In addition, 
the particles exhibited a global vortex fiow. We showed 
that the particle vortices were induced by the convection 
of neutral gas analogously to convective clouds in the 
atmosphere. In turn, the gas convection was triggered 
by the thermal creep along the inhomogeneously heated 
tube surface. 
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